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Abstract 

Event-by-event fluctuations of the kaon to pion number ratio in nucleus-nucleus collisions are studied 
within the statistical hadron-resonance gas model (SM) for different statistical ensembles and in the 
Hadron-String-Dynamics (HSD) transport approach. We find that the HSD model can qualitatively 
reproduce the measured excitation function for the K/n ratio fluctuations in central Au+Au (or 
Pb+Pb) collisions from low SPS up to top RHIC energies. Substantial differences in the HSD and SM 
results are found for the fluctuations and correlations of the kaon and pion numbers. These predictions 
impose a challenge for future experiments. 
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I. INTRODUCTION 



The study of event-by-event fluctuations in high energy nucleus-nucleus (A+A) collisions 
opens new possibilities to investigate the phase transition between hadronic and partonic matter 
as well as the QCD critical point (cf. the reviews By measuring the fluctuations one might 
observe anomalies from the onset of deconfinement 21] and dynamical instabilities when the 
expanding system goes through the 1-st order transition line between the quark-gluon plasma 
and the hadron gas 3]. Furthermore, the QCD critical point may be signaled by a characteristic 
pattern in the fluctuations as pointed out in Ref. ^. However only recently, due to a rapid 
development of experimental techniques, first measurements of the event-by-event fiuctuations 
of particle multiplicities j^, 6, [?, Q] and transverse momenta 9] in nucleus-nucleus collisions 
have been performed. 

From the theoretical side such event-by-event fiuctuations for charged hadron multiplicities 



in nucleus-nucleus colhsions) have been studied in statistical models 



16 



17l . isl . [l9l | and in djTiamical transport approaches 



2Q|, 
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[lO, 
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231], which have been used 



as important tools to investigate high-energy nuclear collisions. We recall that the statistical 



models reproduce the mean multiplicities of the produced hadrons (see e.g. Refs. |24j. |25|. 1261]). 
whereas the transport models (see, e.g., Refs. [271, l28|, l29|]) provide, in addition, a dynamical 
description of the various bulk properties of the system. By studying the various fiuctuations 
within statistical and transport models we have found out that fiuctuations provide an ex- 
tremely sensitive observable - depending on the details of the models - which are partly washed 
out by looking at general quantities such as ensemble averages. 

In particular, there is a qualitative difference in the properties of the mean multiplicity 
and the scaled variance of the multiplicity distribution in statistical models. In the case of 
mean multiplicities the results obtained within the grand canonical ensemble (GCE), canonical 
ensemble (CE), and micro-canonical ensemble (MCE) approach each other in the large vol- 
ume limit. One refers here to the thermodynamical equivalence of the statistical ensembles. 
However, it was recently found 1^ that corresponding results for the scaled variances are 
different in the GCE, CE and MCE ensembles, and thus the scaled variance is sensitive to 
global conservation laws obeyed by a statistical system. These differences are preserved in the 
thermodynamic limit. 

Also there is a qualitative difference in the behavior of the scaled variances of multiplicity 
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distributions in statistical and transport models. The transport models predict 2l|, |22] that the 
scaled variances in central nucleus-nucleus collisions remain close to the corresponding values in 
proton-proton collisions and increase with collision energy in the same way as the corresponding 
multiplicities, whereas in the statistical models the scaled variances approach finite values at 
high collision energy, i.e. become independent of energy. Accordingly, the differences in the 
scaled variance of charged hadrons can be about factor of 10 at the top RHIC energy 21|. Only 
upcoming experimental data can clarify the situation. 

The QGP stage may form a specific set of primordial fiuctuation signals. A well known 
example is the equilibrium electric charge fiuctuation in QGP which is about a factor 2-3 



smaller than in an equilibrium hadron gas 



30, 



31| . To observe primordial QGP fiuctuations 



they should be frozen out during expansion, hadronization, and further hadron-hadron re- 



scatterings. Evolution and survival of the conserved charge fiuctuations in systems 



"ormed in 



nucleus-nucleus collisions at the SPS and RHIC energies were discussed in Refs. [32|, |33|. Note 
that both the statistical models and the HSD approach used in our study do not include the 
quark-gluon degrees of freedom. Thus, the fiuctuations in the QGP are outside of the scope of 
the present paper. 

The measurement of the fiuctuations in the kaon to pion ratio by the NA49 Collaboration js^ 
was the first event-by-event measurement in nucleus-nucleus collisions. It was suggested that 
this ratio might allow to distinguish the enhanced strangeness production attributed to the 
QGP phase. Nowadays, the excitation function for this observable is available in a_wide range 
of energies: from the NA49 collaboration in Pb-|-Pb co 
the STAR collaboration in Au+Au collisions at RHIC 



hsions at the CERN SPS fl and from 
B|. First statistical model estimates of 
the K/t^ fiuctuations have been reported in Refs. 3J, 1351], and results from the transport model 
UrQMD in Ref. [36 1. 

In this paper we present a systematic study of statistical model results (in different ensem- 
bles) in comparison to HSD transport model results for the fluctuations in the kaon to pion 
number ratio. The paper is organized as follows: In Section II the characteristic deflnitions for 
fluctuations in particle number ratios are introduced. In Section HI the relevant formulas of 
the statistical models (in different ensembles) are presented. Statistical and HSD model results 
for the fluctuations in the kaon to pion ratio for central nucleus-nucleus collisions are compared 
in Section IV. In Section V the HSD transport model results are additionally confronted with 
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the available data on K/t{ fluctuations. A summary closes the paper in Section VI. 



II. MEASURES OF PARTICLE RATIO FLUCTUATIONS 
A. Notations and Approximations 

Let us introduce some notations. We deflne the deviation from the average number 

{N a) of the particle species A by = (Na) + AA^^. Then we deflne covariance for species A 
and B 



AiNA,NB) = {ANaANb) = {NaNb) - {Na){Nb) 



scaled variance 



ua 



A {Na, Na) {{ANAf) (iVl) - (Na)' 



and correlation coefficient 



(Na) 



Pab 



(Na) 



{ANa ANb) 



(Na) 



[{{ANAf) {{ANb) 



2\lV2 ■ 



The fluctuations of the ratio Rab = Na/Nb will be characterised by [3J, l35| 



{{^Rab?) 
{RabV 



Using the expansion, 

Na {Na) + ANa {Na) + ANa 



Nb {Nb) + ANb 



{Nb) 



X 



1 



ANf 



+ 



fANB^ ' 



{Nb) \{Nb) 



(2) 



(3) 



(4) 



(5) 



one flnds to second order in ANa/ {Na) and ANb/ {Nb) the average value and the fluctuations 
of the A to B ratio: 



{Rab) 



1 + 



LOb 



AiNA,NB) 

{Nb) {Na){Nb) J 



{Na) 
{Nb) 

A{Na,Na) , AiNB,NB) _ A{Na,Nb) 



(6) 



{Na)' 



UJa ^ ^B 



{Na) {Nb) 



{Nb)' 

2pAB 



UJaUJb 



{Na){Nb) 

1/2 



{Na){Nb) 



(7) 
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If species A and B fluctuate independently according to Poisson distributions (this takes place, 
for example, in the GCE for an ideal Boltzmann gas) one finds ua = = 1 and pab = 0. 
Equation ([7]) then reads 

In a thermal gas, the average multiplicities are proportional to the system volume V. Equation 
(IH]) demonstrates then a simple dependence of cr^ oc 1/V on the system volume. 

A few examples concerning to Eq. ([7]) are appropriate here. When (Nb) ^ {^a), e.g., 
A = + K~ and 5 = vr"*" + vr", the quantity ([7]) is dominated by the fluctuations of less 
abundant particles. When {N a) — (Nb), e.g., A = 7t~^ and B = 7r~ , the correlation term in 
Eq. ([7]) may become especially important. A resonance decaying always into a Tr+vr^-pair does 
not contribute to cr^ ([7]), but contributes to the 7r+ and 7i~ average multiplicities. This leads 
jssi to a suppression of cr^ ([7]) in comparison to its value given by Eq. ([8]). For example, if 
all 7r+ and tt^ particles come in pairs from the decay of resonances, one flnds the correlation 
coefficient p^+jr- = 1 in Eq. (^^, and thus cr^ = 0. In this case, the numbers of vr"'" and vr" 
fiuctuate as the number of resonances, but the ratio tt^/tt^ does not fiuctuate. 

B. Mixed Events Procedure 



The experimental da , 
dynamical fiuctuations 



37| 



or Na/Nb fiuctuations are usually presented in terms of the so called 



adyn = sign((T^ - a^i^) |o-^ - , (9) 

where cr^ is defined by Eq. ([7]) and a.^^^ corresponds to the following mixed events procedure^. 
One takes a large number of nucleus-nucleus collision events and measures the numbers of Na 
and Nb in each event. Then all A and B particles from all events are combined into one set. 
The construction of mixed events is done as follows: One fixes a random number = Na + Nb 
according to the experimental probability distribution P{N), takes randomly A^ particles {A 



^ Other dynamical measures, such as $ [38|,|39| and F [35[, can be also used. 

^ We describe the idealized mixed events procedure appropriate for the model analysis. The real experimen- 
tal mixed events procedure is more complicated and includes experimental uncertainties, such as particle 
identification etc. 
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and/or B) from the whole set, fixes the values of Na and Nb, and returns these particles into 
the set. This is the mixed event number one. Then one constructs event number 2, number 3, 
etc. 

Note that the number of events is much larger than the number of hadrons, N, in any single 
event. Therefore, the probabilities pa and pb = ^ — Pa, to take the A and B species from 
the whole set, can be considered as constant values during the event construction. Another 
consequence of a large number of events is the fact that A and B particles in any constructed 
mixed event belong to different physical events of nucleus-nucleus collisions. Therefore, the 
correlations between the Nb and Na numbers in a physical event are expected to be destroyed 
in a mixed event. This is the main purpose of the mixed events construction. For any function 
f{NA, Nb) the mixed events averaging is then defined as 



(/(iVA,iVB)Ux. = $^P(iV) Yl fi^A^^B) 5{N-NA-NB) ^^;^^_^f p^^p^- . (10) 

N Na,Nb ^' ^' 



The straightforward calculations of mixed averages ( jTOjl can be simplified by introducing the 
generating function Z{x,y), 



Z{x,y) ^ 5^P(iV) J2 m - Na - Nb)^-^^^^^ {xpAf- {ypBf- 

N Na,Nb ^' 

= J2P(N)ixpA + yPBf , (11) 

N 

which depends on auxiliary variables x and y. The averages ffTOl) are then expressed as x- and 
^/-derivatives of Z{x,y) at x = y = 1. One finds: 

{NA)m^. = (1^) = PA {N) , {Nb)^^. = = Pb {N) , (12) 

{NAiNA - 1))™. = (^) = p\ {N{N - 1)) , (13) 

{Nb{Nb-1)U.= ^ =pI{N{N-1)), (14) 

{NANB)mix - {NA)mix{NB)mix = i^g^ J = PaPB {N) , (15) 



where 



(iV) ^ $^iVP(iV), {N') ^ 5^iV^P(iV), UN ^ ^^^^^W^- (16) 

N N ^ ' 
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Calculating the Na/Nb fluctuations for mixed events according to Eq. (JTj) one gets: 



{Na,Na) , A.^i,{NB,NB) 



B 



[{Na) 
1 



+ 



{Na? 
^jv - 1 
{N) 

1 



+ 



{NbY 



+ 



{Nb) 



{N) 



{Na){Nb) 



{N) 



{Na) ^ {Nb) ' ^^^^ 
A comparison of the flnal result in Eq. fll7l) with Eq. (jSj) shows that the mixed events procedure 
gives the same o"^ for Na/Nb fluctuations as in the GCE formulation for an ideal Boltzmann 
gas, i.e. for ua = ^^b = ^ and pab = 0. If = 1 (e.g. for the Poisson distribution P{N)), one 
^mtx _ ]^ g^j^^ pTuix _ Q Otherwise, if 7^ 1, the mixed events procedure 



indeed flnds 00'^^^ 



leads to u;™^ 7^ 1, cu™^ 7^ 1, and to non-zero NaNb correlations, as seen from the second line of 
Eq. ( fT7j) . Thus, if, e.g., event-by-event fluctuations in the total number of pions and kaons are 
stronger than Poissonian ones, i.e. ujn > 1, positive pion-kaon correlations appear in the mixed 
events. They lead to larger (smaller) N^ in the sample of mixed events with larger (smaller) 
A^^. However, the flnal result for a,^^^. ffT7|) is still the same as for cjat = 1, it does not depend 
on the speciflc form of P{N). Non-trivial (cu^J 7^ 1) fluctuations of Na and Nb as well as 
non-zero p^]f correlations may exist in the mixed events procedure, but they are cancelled out 



m cr: 



III. FLUCTUATIONS OF RATIOS IN STATISTICAL MODELS 
A. Quantum Statistics and Resonance Decays 

The occupation numbers, Upj, of single quantum states (with flxed projection of particle 
spin) labelled by the momentum vector p are equal to Upj = 0, 1, . . . , 00 for bosons and 
Up J = 0, 1 for fermions. Their average values are 

exp [{epj - fij) /T\ - aj 

and their fluctuations read 

( (A^pj)^ )gce = {{ripj - {npj)f)g^^ = {upj) (1 + aj (upj)) = - , (19) 

where T is the system temperature, nij is the mass of a particle j, tpj = + m| is the single 
particle energy. The value of aj depends on quantum statistics, i.e. +1 for bosons and —1 for 
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fermions, while aj = gives the Boltzmann approximation. The chemical potential fij of a 
species j equals to: /i-,- = qj /ig + bj + Sj /i^, where the particle electric 

charge, baryon number, and strangeness, respectively, while /ig, /i^, fis are the corresponding 
chemical potentials which regulate the average values of these global conserved charges in the 
GCE. 

In the equilibrium hadron-resonance gas model the mean number of primary particles (or 
resonances) is calculated as: 

p 

where V is the system volume and Qj is the degeneracy factor of a particle of species j (the 
number of spin states). In the thermodynamic limit, V oo, the sum over the momentum 
states can be substituted by a momentum integral. 

It is convenient to introduce a microscopic correlator, (Arip jAn^.i), which in the GCE has 
the simple form: 

{Aupj Ank,i)3ce = j 6ij 5pk . (21) 

Hence there are no correlations between different particle species, i ^ j, and/or between dif- 
ferent momentum states, p 7^ k. Only the Bose enhancement, v^j > (ripj) for aj = 1, and the 
Fermi suppression, v^j < (upj) for aj = —1, exist for fluctuations of primary particles in the 
GCE. The correlator ([1]) can be presented in terms of microscopic correlators fl2Tl) : 

{AN* AN* = Yl <^^P.J- ^^k,.),ce = S,, ""h ■ (22) 

P,k p 

In the case i = j equation (l22l) gives the variance of primordial particles (before resonance 
decays) in the GCE. 

For the hadron resonance gas formed in relativistic A+A collisions the corrections due to 
quantum statistics (Bose enhancement and Fermi suppression) are smalF. For the pion gas 
at T = 160 MeV, one finds cj^ = 1.1, instead of = 1 for Boltzmann particles. The quan- 
tum statistics effects are even smaller for heavier particles like kaons and almost negligible for 
resonances. 



Possible strong Bose effects are discussed in Ref. 12 1 
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The average final (after resonance decays) multiplicities {Ni) are equal to: 

R 



(23) 



In Eq. (l23l) . A^* denotes the number of stable primary hadrons of species i, the summation 
runs over all types of resonances R, and {ni)ji = J2r b^'^fr is the average over resonance decay 
channels. The parameters are the branching ratios of the r-th branches, nf^ is the number 
of particles of species i produced in resonance R decays via a decay mode r. The index r runs 
over all decay channels of a resonance R with the requirement Xlr^r^ ~ ■'■^ GCE the 
correlator ([1]) after resonance decays can be calculated as [35|: 



{ANA/^NB)gce = (AiV* AiV*),,e + [{/\NI) {u a) R^u b) R + (iV«)(An^AnB)fi] , (24) 



R 



where (An^ An^)^ = Y.r^r^A,r^B,r - {n'A)R{nB)R ■ 



B. Global Conservation Laws 



In the MCE, the energy and conserved charges are fixed exactly for each microscopic state of 
the system. This leads to two modifications in comparison with the GCE. First, additional terms 
appear for the primordial microscopic correlators in the MCE. They refiect the ( ant i) correlations 
between different particles, i ^ ?', and different momentum levels, p 7^ k, due to charge and 



energy conservation in the MCE 



1^, 



(Aripj Ank,j)mce = v^j 6ij 5pk TJV^ t ^i^i^ii + hhjMi^i, + SiSjMs 



+ {QiSj + qjSi) Mgs - {qibj + qjhi) Mgb - {hsj + hjSi) M^s 

+ epjCkiMe, - {qiepj + qje^) Mg^ + {bitpj + bje^) Mb^ - {siCpj + Sje^) Mse ] , (25) 
where |A| is the determinant and Mjj are the minors of the following matrix, 

/A(g2) A(6g) A(sg) A(eg)\ 

A(g6) A(62) A(s6) A(e6) 

A(gs) A(6s) A(s2) A(es) 

\A(ge) A(6e) A(se) A(e2)/ 



A 



(26) 
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with the elements, A{q^) = q^v^j , A{qb) = J^pj Qj^j^h > A{qe) = qjepjvlj , etc. 
The sum, , means integration over momentum p, and the summation over all hadron- 

resonance species j contained in the model. The first term in the r.h.s. of Eq. fl25l) corresponds 
to the microscopic correlator (I2T!) in the GCE. Note, that the presence of the terms containing 



the single particle energy epj = y + m'j in Eq. fl25l) is a consequence of energy conservation. 
In the CE, only charges are conserved, thus the terms containing e-pj in Eq. (1251) are absent. 
The matrix A in Eq. ( l26l) then becomes a 3 x 3 matrix (see Ref. 13[). An important property 
of the microscopic correlator method is that the particle number fluctuations and the corre- 
lations in the MCE or CE, although being different from those in the GCE, are expressed by 
quantities calculated within the GCE. The microscopic correlator fl25|) can be used to calculate 
the primordial particle (or resonances) correlator in the MCE (or in the CE): 

{ANi ANj )mce = (^^P.* Ank,j)mce ■ (27) 

p,k 

A second feature of the MCE (or CE) is the modification of the resonance decay contribution 
to the fluctuations in comparison to the GCE fl24|) . In the MCE (or CE) it reads 

{ANAANB)mce = {AN\AN*s)mce + ^^(A^fl) {Aua Aub) R + ^^{ANl AN r) ^ce {us) R 

R R 

+ Y{AN*B ANR)mce Mr + 5^(AAr^ ANr,)^,^ Mr Mr' • (28) 

R R,R' 

Additional terms in Eq. ( l28l) compared to Eq. ( 12^ are due to the correlations (for primordial 
particles) induced by energy and charge conservations in the MCE. Eq. f l28|) has the same form 
in the CE [l^] and MCE [l^, the difference between these two ensembles appears because 
of different microscopic correlators fl25l) . The microscopic correlators of the MCE together 
with Eq. (1271) should be used to calculate the correlators (AA'^AA"^)mce , {AN^ AN r) mce , 
{ANX AN r) mce , {AN% AN r) mce , and (AA^^ ANR,)mcc entering in Eq. The correlators 
(l28ll define finally the scaled variances oja and ujr ^ and correlations pab (E]) between the 
A^ and A"^ numbers. Together with the average multiplicities (A^a) and {Nr) they define 
completely the fluctuations cr^ ([7]) of the AtoB number ratio. 

IV. STATISTICAL AND HSD MODEL RESULTS FOR THE K/t: RATIO 



In this section we present the results of the hadron-resonance gas statistical model (SM) 
and the HSD transport model for the fluctuations of the AT/vr ratio in central nucleus-nucleus 
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collisions. To carry out the SM calculations one has to fix the chemical freeze-out parameters. 
The dependence of the chemical potential /i^ on the collision energy is parameterized as I25I: 
fJ'B (a/stvtv) = 1.308 GeV • (1 + 0.273 y^s/vw)^^ , where the center of mass nucleon-nucleon 
collision energy, ^/snn, is taken in units of GeV. The system is assumed to be net strangeness 
free, i.e. S = 0, and to have the charge to baryon ratio of the initial colliding nuclei, i.e. 
Q/B = 0.4. These two conditions define the system strange, fis, and electric, fig, chemical 
potentials. For the chemical freeze-out condition we chose the average energy per particle, 
(E) I (N) = 1 GeV 40]. Finally, the strangeness saturation factor, 75, is parameterized as in 
Ref. 26|: 75 = 1 — 0.396 exp (— 1.23 T/^b) ■ This determines all parameters of the model. 



An extended version of the THERMITS framework 1:41.] is used for the SM calculations (for 
more details see Ref. |l4]). Note that average multiplicities of pions and kaons measured in 
central nucleus-nucleus collisions at SPS and RHIC energies 42| are nicely described in the SM 
(see, e.g., Refs. ybej) as well as HSD (cf. Ref. |29|). 



A. Results for uj^^ ^ni and pxn 



2 , ^TT r, 



(29) 



According to Eq. ([7]) the fiuctuations of the K = + K to vr = 7r+ + vr ratio is given by 

1 1/2 

The values of uj^^, ujk and pKn in different statistical ensembles are presented in Table I and 
for the HSD simulations of Pb+Pb (Au+Au) central (with impact parameter 6 = 0) collisions 
in Table II. Both the SM and HSD results are shown in Figs. [Hand [H Let us first comment the 
SM results. In the SM the scaled variances u.,, and ojk and correlation parameter pk-k approach 
finite values in the thermodynamic limit of large volumes. These limiting values are presented 
in Table I and in Fig. [1] and [2l For central Pb+Pb and Au+Au collisions the corresponding 
volumes in the SM are large enough. Finite volume corrections are expected to be on the level 
of a few percent. The finite volume effects for the scaled variances and correlation parameters 
in the CE and MCE are, however, difficult to calculate (see Ref. 15|) and they will not be 



considered in the present paper. The GCE values of cj^ and ujk reflect the Bose statistics of 
pions and kaons and the contributions from resonance decays. 

The TT-K correlations px-w are due to resonances having simultaneously K and vr mesons in 
their decay products. In the hadron-resonance gas within the GCE ensemble, these quantum 
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TABLE I: The chemical freeze-out parameters T and /ib for central Pb+Pb (Au+Au) colUsions at 
different cm. energies ^/SNN■ The hadron-resonance gas model results are presented for final (after 
resonance decays) number densities of pions n^^ and kaons nx (they are the same in all statistical 
ensembles), scaled variances Wtt, ojk and correlation parameter px-w in the GCE, CE, and MCE. 
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TABLE IL The HSD average multiphcities (A^tt), {Nk) and values of lOt^, ^k, and pKn for central 
(impact parameter 6 = 0) Pb-|-Pb (Au-|-Au) collisions at different cm. energies ^/sNN. 
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T 1 1 — I — I 1 1 1 1 1 1 1 — I — I r 




10 , 100 

FIG. 1: (Color online) The SM results in the GCE, CE, and MCE ensembles and the HSD results 
(impact parameter 6 = 0) are presented for the scaled variances uj^^, ujk for Pb+Pb (Au+Au) collisions 
at different cm. energies ^/sJiN- For comparison the HSD results for inelastic proton-proton collisions 
are also presented in terms of the dotted lines with open circles. 
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FIG. 2: (Color online) The SM results in the GCE, CE, and MCE ensembles and the HSD results 
(impact parameter 6 = 0) are presented for the correlation parameter pKn for Pb+Pb (Au+Au) 
collisions at different cm. energies y^J/vJV- For comparison the HSD results for inelastic proton-proton 
collisions are also presented by the dotted line with open circles. 

statistics and resonance decay effects are responsible for deviations of ujk and from 1, and 
of pktt from 0. The most important effect of an exact charge conservation in the CE ensemble 
is a suppression of the kaon number fluctuation. This happens mainly due to exact strangeness 
conservation and is reflected in smaller CE values of ux at low collision energies in comparison 
to those from the GCE ensemble. The MCE values of ujk and cUtt are further suppressed in 
comparison those from the CE ensemble, which is due to exact energy conservation. The effect 
is stronger for pions than for kaons since pions carry a larger part of the total energy. An 
important feature of the MCE is the anticorrelation between A^^ and N^, i.e. negative values 
of pktt- This is also a consequence of energy conservation for each microscopic state of the 
system in the MCE [l^. The presented results demonstrate that global conservation laws are 
rather important for the values of tu^, ujk, and pK-w- In particular, the exact energy conservation 
strongly suppresses the fluctuations in the pion and kaon numbers and leads to uk < 1 and 
< 1 in the MCE ensemble instead of ujk > 1 and a;^ > 1 in the GCE and CE ensembles. 
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The exact energy conservation changes also the rc-K correlation into an anticorrelation: instead 
of pktt > in the GCE and CE ensembles one finds pktt < in the MCE. 

As seen from Figs. [1] and [2] the HSD results for uj^,, ujk, and pKn (solid lines) are rather 
different from those in the SM. For a comparison the HSD results for inelastic proton-proton 
collisions are also presented in Figs. [T]and[2] (dotted lines). The HSD scaled variances and 
cuk increase at higher energies. A similar behavior has been observed earlier in Ref. [21| for 
the scaled variance of all charged hadrons. The HSD calculations reveal the anticorrelation 
between A^^^ and N^, i.e. negative values of pxn, for low SPS energies, where the infiuence of 
conservation laws is more stringent. 

Comparing this result with the SM (in different ensembles) one may conclude that negative 
values of pKir in HSD appear because of a dominant role of energy conservation in joint tt- 
K production at small collision energies. The HSD values of pktt become, however, positive 
and strongly increases with increasing collision energy. This is due to the contribution of 
heavy strings to joint n-K (or K*) production at high energies in the HSD simulations. Note 
that the HSD results for Uj,, lj^, and pxn in nucleus- nucleus collisions become larger than 
those in proton-proton inelastic reactions at high collision energies. This is due to an increase 
of secondary (i.e. meson-baryon and meson-meson) collisions at higher bombarding energy. 
Thus, a strong deviation of HSD from the SM with increasing energies is a consequence of non- 
equilibrium dynamics in the hadron-string model which is driven by the formation of heavy 
strings and their decay. Indeed, future experimental data on the fluctuations of K, n and Ktt- 
correlations will allow to shed more light on the equilibration pattern achieved in heavy-ion 
collisions at RHIC energies. 

Two comments are appropriate here. The flrst one concerns a correspondence between the 
HSD and SM results. In HSD three charges - net baryon number B (equal to the number of 
participating nucleons), net electric charge Q (equal to the number of participating protons), 
and net strangeness S (equal to zero) - are conserved exactly during the system evolution. 
However, B and Q can fluctuate from event to event because of the fluctuations in the number 
of nucleon participants. They also cause fluctuations of the energy of produced hadrons in the 
HSD simulations. Besides, an essential part of the system energy is transformed to collective 
motion. Thus, even in the sample of the HSD events with flxed number of participants, the 
thermal energy of the created particles can fluctuate from event to event. Both the charge and 
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energy fluctuations in HSD are not of tliermal origin. Tfierefore, an attempt to interpret tlie 
HSD multiplicity fluctuations in statistical terms would require to use a more general concept 
of statistical ensembles with fluctuating extensive quantities jl?!. If^. In particular, large values 
of the scaled variances, uJi ~ (A^j), in high energy proton-proton collisions are also present in 
the HSD simulations of nucleus-nucleus collisions. In the SM model this would require a special 
form of scaling volume fluctuations as recently suggested in Ref. [3]. 

Our second comment concerns the physical origin of the correlation parameter pk-k- Two 
sources of the n-K correlations are: resonance, string decays and electric charge conservation. 
To estimate their relative weights, one can beneflt from measuring the correlations pk-k in the 
separate charge channels: it~K~ and 'k~K~^ as suggested in Ref. 43|. The resonances decaying 
into n~K~^ produce the corresponding correlation, while an analogous correlation in the 7i~K~ 
system is absent. Note that electric charge conservation leads also to qualitatively different 
correlation effects in 7c~K~ and tt~K~^ channels. 



B. Results for a, amix, and adyn 

The fluctuation in the kaon to pion ratio is dominated by the fluctuations of kaons alone since 
the average multiplicity of kaons is about 10 times smaller than that of pions. Thus, the 1-st 
term in the r.h.s. of Eq. (!29l) gives the dominant contribution, while the 2-nd and 3-rd terms 
in fl29|) give only small corrections. The model calculations of fl29l) require, in addition to ujk, 
cUtt, and pktt values, the knowledge of the average multiplicities {Nk) and {Nt,). For the HSD 
simulations (impact parameter 6 = in Pb-|-Pb collisions at SPS energies and Au-|-Au collisions 
at RHIC) the corresponding average multiplicities are presented in Table 11. To flx average 
multiplicities in the SM one needs to choose the system volume. For each collision energy we 
fix the volume of the statistical system to obtain the same kaon average multiplicity in the SM 
as in the HSD calculations: {NK)stat = {Nk)hsd- We recall that average multiplicities of kaons 
and pions are the same in all statistical ensembles. The SM volume in central Pb-|-Pb (Au-|-Au) 
collisions is large enough and all statistical ensembles are thermodynamically equivalent for the 
average pion and kaon multiplicities since these multiplicities are much larger than 1. 

In Fig. [3] the values of cr (in percent) - calculated according to Eq. fl29|) and Eq. f|T7|) - are 
presented in the left and right panel, respectively, for the SM in different ensembles as well 
as for the HSD simulations. The first conclusion from Fig. [3] (left) is that all results for a in 
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the different models are rather similar. One observes a monotonia decrease of a with collision 
energy. This is just because of an increase of the kaon and pion average multiplicities with 
collision energy. The mixed event fluctuations (ymix in the model analysis are fully defined by 
these average multiplicities according to Eq. f|T7j) . The values of (Jmix are therefore the same 
in the different statistical ensembles. They are also very close to the HSD values because we 
have fixed the statistical system volume to obtain the same kaon average multiplicities in the 
statistical model as in HSD at each collision energy. As seen from Yig.^{right) the requirement 
of {NK)stat = {Nk)hsd leads to practically equal values of a^ix in both HSD and the SM. 
18 
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FIG. 3: (Color online) Left: The SM results in the GCE, CE, and MCE ensembles as well as the HSD 
results (impact parameter 6 = 0) are presented for a-100% defined by Eq. ([29]) for Pb+Pb (Au+Au) 
collisions at different cm. energies ^/sNN■ Right: The same as in the left panel, but for amix-WOVo in 
mixed events defined by Eq. ([IT]), a^j^. = 1/{Nk) + 1/{N^). 

Differences between the statistical ensembles as well as between the statistical and HSD 
results become visible for other measures of K/tt fluctuations such as adyn deflned by Eq. (Q 
and F = cr^/cr^j^. They are shown in Fig. H] left and right, respectively. At small collision 
energies the CE and MCE results in Fig. H] demonstrate negative values of <7dyn, respectively 
F < 1. When the collision energy increases, adyn in the CE and MCE ensembles becomes 
positive, i.e. F > 1. Moreover, the different statistical ensembles approach to the same values 
of (Jdyn and F at high collision energy. In the SM the values of a and amix approach zero at high 
collision energies due to an increase of the average multiplicities. The same limit should be also 
valid for a^yn in the SM. In contrast, the measure F shows a different behavior at high energies: 
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FIG. 4: (Color online) Left: The results for the K/n fluctuations at different c .m. energies -y/sTViv 
the GCE, CE, and MCE ensembles as well as from HSD (impact parameter 6 = 0) are presented for 
<7d2yn"100% defined by Eq. Q. Right: The same as in the left panel but for F = cr^/cr^^^. 

the SM gives F = 1.05 in the high energy limit, while the HSD result for F demonstrates a 
monotonic increase with collision energy. An interesting feature of the SM is approximately the 
same result for a (and, thus, adyn and F) in the CE and MCE ensembles. From Table I and 
Figs. [Hand [2] one observes that both ujk, ^-k and px-w are rather different in the CE and MCE. 
Thus, as discussed above, an exact energy conservation influences the particle scaled variances 
and correlations. These changes are, however, cancelled out in the fluctuations of the kaon to 
pion ratio. 



C. Volume Fluctuations 



It has been mentioned in the literature (see, e.g., Ref. [3J|) that the particle number ratio 
is independent of volume fluctuations since both multiplicities are proportional to the volume. 
In fact, the average multiplicities (Nk) and (A'^^), but not and N^^, are proportional to 
the system volume. Let us consider the problem in the SM assuming the presence of volume 
fluctuations at fixed values of T and fiB- This assumption corresponds approximately to volume 
fluctuations in nucleus-nucleus collisions from different impact parameters in each collision 
event. Under these assumptions the SM values in Table I remain the same for any volume (if 
only this volume is large enough and the finite size corrections can be neglected). However, the 
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average hadron multiplicities are proportional to the volume. Therefore, the SM result for o"^ 
reads, = ctqVo/V, where Vq is the average system volume, and ctq is calculated for the average 
multiplicities corresponding to this average volume Vq. Expanding Vq/V = Vq/{Vo + 6V) in 
powers of 6V/V0, one finds to second order in 6V/V0, 



^2 rNJ 2 



1^2 



(30) 



where 

{{dVf) = j dV {V - Vof W{V) (31) 

corresponds to an average over the volume distribution function Vr(l^) which describes the 
volume fluctuations. As clearly seen from Eq. fl5U]l the volume fluctuations influence, of course, 
the K/tc particle number fluctuations and make them larger. Comparing the K/tt particle 
number fluctuations in, e.g., 1% of most central nucleus-nucleus collisions with those in, e.g., 
10% one should take into account two effects. First, in the 10% sample the average volume 
Vq is smaller than that in 1% sample and, thus, Uq in Eq. (l30ll is larger. Second, the volume 
fluctuations ( 13T1) in the 10% sample is larger, and this gives an additional contribution to 
according to Eq. fl30|) . 

One may also consider volume fluctuations at fixed energy and conserved charges (see, e.g., 
Ref. [3]). In this case the connection between the average multiplicity and the volume becomes 
more complicated. The volume fluctuation within the MCE ensemble can strongly affect the 
fluctuations in the particle number ratios. This possibility will be discussed in more detail in a 
forthcoming study. 



V. EXCITATION FUNCTION FOR THE K/tt RATIO: COMPARISON WITH DATA 



A comparison of the SM results for K/tt fluctuations in different ensembles with the data 
looks problematic at present. This is because of difficulties with implementing the experimental 
acceptance in the SM (see a discussion of this point in Ref. 16Q. A similar problems exist in the 
SM with chemical non-equilibrium effects discussed in Ref. 4^. The experimental acceptance 
can be taken into account in the transport code. In order to compare the HSD calculations 
with the measured data the experimental cuts are applied for the simulated set of the HSD 
events. In Fig. [5] the HSD results for the excitation function in adyn & for the K/ n ratio is 
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shown in comparison with the experimental data measured by the NA49 Collaboration at the 



SPS CERN 



7| and by the STAR Collaboration at BNL RHIC 



For the SPS energies we used a cut piab > 3 GeV/c applied by NA49 to provide a precise 
parfcle identfflcation. For the RHIC e.erg,es the cuts are i,r pseudorapidity. |,| < 1. and ,„ 

the transverse momentum, 0.2 < px < 0.6 GeV/c, [8]. We note also, that the HSD results 
presented in Fig. [5] correspond to the specific centrality selections as in the experiment - the 
NA49 data correspond to the 3.5% most central collisions selected via the veto calorimeter, 
whereas in the STAR experiment the 5% most central events with the highest multiplicities in 
the pseudorapidity range \ri\ < 0.5 have been selected. 

The HSD results - within the acceptance cuts - are shown in Fig. [5] {left) as solid lines. In 
addition in the l.h.s. of Fig. [5] the result for the full acceptance is indicated by a dotted line. 
One can see that the experimental cuts lead to a systematic increase of adyn, however, do not 
change the shape of the excitation function. By comparing the full acceptance line to those in 
Fig. m (left) for 6 = one sees also a small enhancement of adyn which is due to slight decrease 
of the hadron multiplicities and, correspondingly, increase of those fiuctuations. 

I,r the n,Mpa,>e. of F.g.^the HSD results for w.th.u the experimeuta. acceptauce are 
compared with two different versions of UrQMD vl.3 simulations |36l . I45l | and the NA49 data 
in the SPS energy range. The remaining differences between the UrQMD vl.3 calculations 
from 2006 and 2004 at 160 A GeV can be attributed to the differences in implementation of 
acceptance cuts (cf. discussion in Ref. |45]). One sees that the UrQMD model gives practically a 
constant (Xdyn, which is about 40% smaller than the results from HSD at the lowest SPS energy. 
This difference between the two transport models might be attributed to different realizations 
of the string and resonance dynamics in HSD and UrQMD: in UrQMD the strings decay first 
to heavy baryonic and mesonic resonances which only later on decay to 'light' hadrons such 
as kaons and pions. In HSD the strings dominantly decay directly to 'light' hadrons (from 
the pseudoscalar meson octet) or the vector mesons p, u and K* (or the baryon octet and 
decouplet in case of baryon number ±1). As discussed in the previous section, cxdyn is indeed 
very sensitive to the model details at low bombarding energies: the SM in different ensembles 
and the HSD give rather different behavior at the low SPS energies (cf. Fig. HI left). 

While the UrQMD results are available presently only up to the top SPS energy, the HSD 
model shows a good agreement with the recent STAR data [8|] (cf. Fig. [5], left). A good 
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FIG. 5: (Color online) Left: The HSD results for the excitation function in (Jdyn & for the K/tt ratio 
for full acceptance ( dotted line) and within the experimental acceptance (solid line) in comparison to 
the experimental data measured by the NA49 Collaboration at the SPS CERN [7| and by the STAR 
Collaboration at BNL RHIC [8]. 3.5% most central HSD events were selected for the analysis for SPS 
energies and 5% - for the RHIC energies. Right: The HSD results (circles) and two different versions 



of UrQMD (triangles) calculations 



36 



45l | for adyn versus the NA49 data. Statistical uncertainties in 



the transport calculations are shown by error bars. 



agreement with the STAR data [46| for K/tt ratio fluctuations in Cu+Cu at ^/snn =200 GeV 
was also obtained in the Multi-Phase Transport Model (AMPT) [47^ . This is in contrast to the 
corresponding result from the Heavy-Ion- Jet-Interaction Generator (HIJING) model 



over- pre dicts substantially the experimental data 46|]. The difference has been attributed in 



4M which 



Ref. 



46| to an absence of the flnal re-scattering in HIJING which is incorporated in AMPT as 



well as in HSD. 



VI. SUMMARY AND CONCLUSIONS 



We have studied the event-by-event fluctuations of the kaon to pion number ratio in central 
Au+Au (or Pb+Pb) collisions from low SPS up to top RHIC energies within the statistical 
hadron-resonance gas model for different statistical ensembles - the grand canonical ensemble 
(GCE), canonical ensemble (CE), and micro-canonical ensemble (MCE) - and in the Hadron- 
String-Dynamics transport approach. We have obtained substantial differences in the HSD and 
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statistical model results for the scaled variances ujk, and the correlation parameter pKir as 
presented in Figs. [Hand El Thus, the second moments of the multiplicity distributions may serve 
as a good probe for the amount of equilibration achieved in central nucleus-nucleus collisions. 



Note that the differences between the transport and statistical mode^ 
fluctuations and correlations increase with collision energy (see Refs. 



results for multiplicity 
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22|). There are also 



arguments that the behavior of higher moments of event-by-event multiplicities may serve as 



an important signature of the QCD critical point [49l|). 

The observable (Jdyni which characterizes the fluctuations of the kaon to pion ratio, shows 
to be rather sensitive to the details of the model at low collision energies. The CE and MCE 
results in Fig. H] demonstrate negative values for CTdyni while the GCE gives approximately a 
constant positive value for adyn- The HSD results correspond to larger values of adyn than those 
in the GCE statistical model. They even show an increase at lower SPS energies. When the 
collision energy increases, the quantity adyn in the CE and MCE becomes positive. Moreover, 
the different statistical ensembles approach to the same values of (Tdyn at high collision energy. 
This is just because the values of a and (Jmix approach zero at high collision energies. Thus, 
the same limit equal to zero should be also valid for adyn in the statistical models. On the 
other hand, the measure F = a^/a'^j^^ shows another behavior at high energies. The statistical 
models give a constant value F = 1.05 in the high energy limit, while the HSD results for F 
demonstrate a monotonic increase with collision energy. 

We find that the HSD model can qualitatively reproduce the measured excitation function 
for the K/tt ratio fluctuations in central Au+Au (or Pb+Pb) collisions from low SPS up to top 
RHIC energies. We have shown that accounting for the experimental acceptance as well as the 
centrality selection has a relatively small influence on adyn and does not change the shape of 
the adyn excitation function. We conclude, that the HSD hadron-string model - which does not 
have a QGP phase transition and not explicitly includes the quark and gluon degrees of freedom 
- can reproduce qualitatively the experimental excitation function. In particular, it gives the 
rise of adyn with decreasing bombarding energy. This fact brings us to the conclusion that the 
observable enhancement of adyn at low SPS energies might dominantly signal non-equilibrium 
string dynamics rather than a phase transition of hadronic to partonic matter or the QCD 
critical point. 
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